The radiative decay of tp{3770) into the predicted scalar state X(3700) 
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We calculate the radiative decay width of the ^/)(3770) into the dynamically generated scalar 
resonance X(3700) which is predicted in a previous paper. The results show that it is possible that 
the upgraded BEPC-II facility will generate enough statistics in order to observe this decay and 
thus confirm the existence of the X(3700). 
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I. INTRODUCTION 



The extension of SU{3) Lagrangians to SU{4:), with 
the appropriate SU (4) flavor symmetry breaking taken 
into account, has been proved useful in order to describe 
the structure of many charmed resonances as dynamically 
generated states [U, [l] ■ Although it presents a plausible 
explanation for many observed states, its spectrum shows 
many new states which have not been observed. The sta- 
bility of some results, and the good fit achieved between 
the models and the observed experimental states, make 
us confident that the stable predicted states should exist. 
One should be careful with other states that seem to be 
more model dependent. For instance in [1, 0, [E, @| two 
open charm narrow sextets, an axial and a scalar, are 
predicted while in our works [J, 01 , which take into ac- 
count the differences between the meson decay constants 
of charmed and light mesons, these two sextets become 
very broad, the scalar one probably too broad to have 
any experimental relevance. 

The model of also predicts a narrow hidden charm 
scalar resonance around 3.7 GeV and two axial vector res- 
onances, nearly degenerated in mass, one with C-parity 
positive which corresponds to the X{3872), and another 
one with negative C-parity. In order to stimulate the ex- 
perimental search for these states, in the present work 
we present a calculation of the radiative decay width of 
the V'(3770) into the predicted scalar state around 3.7 
GeV. We also discuss the possibility to observe this de- 
cay at the planned BEPC-II facility of Beijing based on 
the expected rate of i/'(3770) production at this facility. 

The work is organized as follows: In the next section 
we present the model for generating scalar resonances 
dynamically. We have done some improvements with re- 
spect to the work of [l|, that we explain. Section III 
is dedicated to the framework to calculate the radiative 
decay width of vector mesons into scalar mesons and in 
section IV we point out our final remarks. 
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II. MODEL 

First we build a field for the mesons belonging to the 
15-plet of SU{4:), with help of the SU{4:) generators, A^: 
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for this field a current is build: 



Ju 



(2) 



note that adding a diagonal matrix (a SU (4) singlet field) 
to $ would give null contribution to J^. In our previous 
work we have used only the mathematical states 778 and 
ryi5 that we have identified with the physical states 77 
and rjc respectively. In the present work we are going to 
consider the mixing between these states and a singlet 
and use the matrix $ in the physical basis: 
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The Lagrangian we use is the current in eq. ([2]) coupled 
to itself plus a mass term: 
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For the mass term we take the matrix M 



diagonal {'m^,rn^,2mj^ — m^,2m|, — m^). This La- 
grangian is a straightforward generalization to SU (4) of 
the usual lowest order chiral Lagrangian used for SU (3) 
in several works 0, Q • Since SU (4) is not a good symme- 
try in nature we are going to explicitly break it in three 
ways. First, since the mass term is not proportional to 
the identity matrix, it already has some amount of SU (4) 
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symmetry breaking. The parameter / in Q is the me- 
son decay constant, we are going to use / = /jr = 93 
MeV for light mesons, and / — = 165 MeV for heavy 
ones, and last we are going to suppress terms in the La- 
grangian where the interaction is driven by the exchange 
of a heavy vector meson, in the vector meson dominance 
picture. For details in this suppression one should refer 
toP,i,i,[l3. 

From the resulting Lagrangian, using usual Feynman 
rules, the transition amplitudes between any possible two 
meson initial and final states that span a coupled channel 
space are calculated, projected in s-wave, and collected 
into a matrix, V. This matrix is then used as the kernel 
to solve the Bethe-Salpeter equation for the T-matrix, 
that in the on-shell formalism of [HI, assumes an 
algebraic form: 

T = V + VGT, (4) 

in this equation G is a diagonal matrix with each one of 
its elements given by the loop function for each channel 
in the coupled channel space. 

The imaginary part of the loop function ensures that 
the T-matrix is unitary, and since this imaginary part is 
known, it is possible to do an analytic continuation for go- 
ing from the first Riemann sheet to the second one. Pos- 
sible physical states (resonances) are identified as poles 
in the T-matrix calculated in the second Riemann sheet 
for the channels which have the threshold below the res- 
onance mass. 

Considering all the mesons in the 15-plet of SU (4) plus 
the singlet field, it is possible to couple 15 channels to 
the quantum numbers C=0 and S=0: tt+tt^, K^K^, 

TT°72^°r]', D + D-, TJcTT°, TT°TT°, K°K°, TJIJ, T^T]' , r( r( , 

D^D'^^ DfDj, ricTj and r/c?y'. Three of these channels 
are pure isospin 1 channels and, since as discussed in pj 
the hidden-charm scalar resonance is an isospin reso- 
nance, it will not couple to those channels. 

TABLE I: Couplings of the pole at (3722-il8) MeV to the 
channels. 



Channel 


Re{gx) [MeV] 


Im{gx) [MeV] 


\gx\ [MeV] 
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22 


VV 


1023 


242 


1051 




1680 


368 


1720 




922 


-417 


1012 


D'^DO 


5962 


1695 


6198 


DTD7 


5901 


-869 


5965 


VcV 


518 


659 


838 




405 


9 


405 



Solving the coupled channels problem for calculating 
the T-matrix with the same parameters as in [l| , but with 



the physical basis for the interacting fields (considering 
the 77, 77' mixing and the rjc as a pure cc state) we find 
now the X(3700) pole at =(3722-il8) MeV. In Table 
H] we show the couplings of this resonance to all channels 
calculated via the residues of the pole to the channels. 

Note that now the X(3700) has a much bigger width 
than in [l| where this width was smaller than 1 MeV. 
This is due to the fact that in [1] the only decay chan- 
nels which were not strongly suppressed by the dynam- 
ics of the interaction were channels involving only heavy 
mesons. The lighter channels were suppressed by the in- 
teraction and the channel 7715778, which is closer in mass to 
the resonance was suppressed because it is a pure SU (3) 
octet while the X is a 5/7(3) singlet. Now, the t] and 77' 
have an important SU (3) singlet component and so the 
couplings of the resonance to channels involving these 
mesons is bigger and therefore the resonance acquires 
a bigger width, although still small compared to other 
hadronic resonances. 

In the next section we present the formalism to evalu- 
ate the radiative decay width of a vector into a scalar. 



III. RADIATIVE DECAY 

We are going to consider the following radiative decay: 

7A(P,e^)->X(Q)-f 7(ir,e^) (5) 

where ip is the second radial excitation of the charmo- 
nium spectrum, the ■0(3770) and X is the dynamically 
generated resonance X(3700). 

We are following here the same steps as in [13] for the 
calculation of the radiative decay width, which is based 
in the formalism used in [ij, [l^ ■ In the reaction ([5]) one 
has two independent four momenta since P = Q + K. 
We chose to work with P and K. Moreover both vector 
particles fulfill the Lorentz condition: 

= (6) 
= (7) 

The amplitude for the decay in ([5]) will be given by 

-tM = ey^T,,, (8) 

and since the problem has two independent four mo- 
menta, by Lorentz invariance one may write 

V = ag^^ + bP^P^ + cP^K^ + dP^K^ + eK^K^{9) 

Now one can realize that by Lorentz condition (H)),!!!]) the 
terms with coefficients b, c and e will not contribute to 
the radiative decay amplitude. 
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FIG. 1: Diagram that contains the d term. 
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^0 J (Stt) 
{q' + yK)^{q' + il-x)P),, 



(16) 
(17) 



(18) 



where in the last equation we have made a change of 
variables q = q' + {1 — x)P + yK and we have defined s: 



(1 - x){xrni - ml- 2yP.K) 



(19) 



Now we can identify the d term in the expression of 
equation and perform the q' integral. 



Applying the gauge condition [K^Tf^,, = 0) to the ex- 
pression in ^ one obtains that the two remaining coeffi- 
cients are related: a = —dK.P, so one needs to calculate 
only one of them in order to obtain the full gauge invari- 
ant amplitude for the process. We are going to calculate 
the d term, which comes from only one diagram, illus- 
trated in figure[TJ In the figure mi and TO2 are the masses 
of the charged mesons to which the X{3700) can couple, 
we will consider only the D^D^ and DfD^ channels, 
since the other charged channels to which it couples have 
negligible couplings compared to these ones, see table HI 

Given the basic couplings of ip{P) {q)D^ [P — q): 



of X(3700) to PP: 



-igi,e^Jq-P + q)^, (10) 



-iMx^pp = -igx, 
and of the photon to PP: 



iM 



-f^pp 



-ieeUq- 



(11) 



(12) 



for a positive pseusdoscalar meson and e being the mod- 
ulus of the electric charge, the amplitude for this process 
is given by: 



iM 



d'^q 
(2^ 
1 



[P — qY — m^ (f- — m\ [q — K)^ — 
Using Feynman parameterization one gets: 



5(13) 



1 

abc 

a = (P-q) 



dx / dy ^ 

Jo (a + {b — a)x + {c — b)yy 



(14) 
(15) 



/ 



d'^q' 
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■ief 



2 s + ie 



(20) 



The coupling of the photon to the other particle (to2 
in figure [T]) gives, in the present case, an identical contri- 
bution since D~ have opposite signs and equal mass 
and the same occurs for DfD~. Hence for each of these 
two channels we obtain: 



. gj^gxe 
' 2^2 



dx 

"'0 



dy 



s -\- ie 



(21) 



and we must sum this expression for the two channels 
D^D^ and DfD^ using the coupling g^ to each channel 
that we evaluate below, and gx which appears in Table 

in 

With the d term one can calculate the a term and 
with these two terms the full gauge invariant amplitude 
A4 for the radiative decay. Making an average over the 
possible polarizations of the 1/^(3770) state of the absolute 
value square of this amplitude, summing over the possible 
photon polarizations and integrating over phase-space, 
one obtains the radiative decay width: 



P 



iPKf\d\', 



(22) 



where p is the relative center of mass three momentum 
of the two particles in the final state. 

There is still one parameter to be calculated which is 
the coupling of the vector -0(3770) to the two possible 
channels. 

At tree level the decay width of a vector meson into 
two pseudoscalars is given by: 



1 



V-^PP 



2„3 



87rM2 



(23) 



where g is the VPP coupling and p is the center of mass 
relative momentum of the two pseudoscalars in the final 
state. 
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The total decay width of the ^lJ{3770) is 25 MeV where 
36% of this width is coming from D^D^ final state, ac- 
cording to [iBl ■ Using 9 MeV for F in eq. ([23]) one obtains 
11.7 for the value of g. The ^^(3770) is the 35* radial ex- 
citation of the J/ip vector meson. Hence we expect the 
flavor dynamics of this state to be identical to that of 
the J/ip. This means that the couplings of -0(3770) to 
the different PP states will be the same as those of the 
J/f/; up to a global proportionality constant related to 
the overlap of the radial wave functions. 

We follow the work of [l| and write the VPP La- 
grangian in St/ (4) formalism: 



Cppv = 



V2 



rr([a^$,$]i/^) 



(24) 



with 



V2 



p 

V2 



K 
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K*- 
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_ A'*" D* 

V2__V2 



D*+ 



d: 



D*- 
J/V' / 



This guarantees that the coupling of -0(3770) to D^D^ 
and to DfDj, which are needed in the loop of figure [U 
are equal, like in the case of J/ip- 

With the expressions in eqs. (PT|) and ([^ it is 
straightforward to evaluate the radiative decay width of 
the ^'(3770) into the X(3700) resonance. 

For the masses of the mesons we have used: m£)=1868 
MeV, mD^=1965 MeV, tox=3722 MeV and m^=3772 
MeV. 

Applying the expressions in (PT|) and one obtains 
for the radiative decay width of the -0(3770) into the 
X(3700): 



0.65 KeV. 



(25) 



Taking into account the finite width of the -0(3770) or 
of the X(3700) in the evaluation of this radiative decay 
has less than 10% effect in the result. 

The value in (^5]) corresponds to a branching fraction 

of: 



B(V'(3770) ^ X(3700) + 7) 
B(V'(3770) anything) 



= 2.6 X 10" 



(26) 



The branching ratio is of the order of magnitude of the 
ao (980)7 or /o (980)7 decays [ll [Hi, which 
proceed with similar loops but involving kaons instead of 
D mesons 19.,^20,]. 

It is interesting to make an estimative of the theoreti- 
cal uncertainties. In [2l| it was suggested that the bump 
in the DD spectrum close to threshold observed at Belle 
in the e+e" J/^DD could be due to the exis- 
tence of the hidden charm state that we have considered. 



The alternative explanation as a new state above the DD 
threshold, mildly suggested in , had a value larger 
than the one obtained assuming the peak as due to the 
hidden charm state below threshold. The best fit to the 
data was obtained with a mass of the X{370Q) around 
3720 MeV. We have changed the subtraction constant of 
the loop function in order to generate the X(3700) pole 
at different positions which still produced good values 
in the fit of the e+e" J/tfiDD data, and we get the 
results for the radiative decay width shown in Table [III 



TABLE II: Results for different pole positions of the X(3700) 
resonance. The value in the E-, column corresponds to the 
photon momentum. 



Oh 


VSr,oie [MeV] 


go [GeV] 


go, [GeV] 


E-, [MeV] 


r [KeV] 


-1.40 


3702-127 


8.064-11.29 


7.83-11.34 


96.73 


3.24 


-1.35 


3713-i21 


6.524-11.54 


6.39-11.06 


58.54 


0.99 


-1.30 


3722-118 


5.964-11.69 


5.90-10.87 


49.67 


0.65 


-1.25 


3730-115 


5.39-H11.91 


5.42-10.60 


41.77 


0.43 



The feasibility of the experiment at BES relies u pon 
the statistics of ?/'(3770) production. According to [2^, 
m in BEPC-II it is expected to have 3.8 x 10^ 0(3770) 
events in one year of run, which would correspond to 
around 1000 events into V(3700)7 for F^^X7=0.65 KeV, 
with the photon energy peaking around 50 MeV. For 
other values of the V(3700) mass the photon peak po- 
sition changes and so does the rate, but in all cases dis- 
regarding technical problems that are beyond our reach, 
the statistics of one year is far more than sufficient to 
determine this peak with the required precision. 

To further support our claim we would like to quote 
that CLEO 25] has reported the observation of the 
-0(3770) decay into 7Xco(l^) which has a decay width 
of 172 KeV, about two orders of magnitude bigger than 
our resuhs for 7^(3700). Since BEPC-II will have two 
orders of magnitude more statistics than CLEO, the mea- 
surement that we suggest seams to be at reach. 



IV. OVERVIEW 

The fact that the predicted hidden charm resonance 
X(3700) lies below the ?/'(3770), together with the plans 
to produce this latter particle at BEPC-II, opens up the 
possibility to see the V(3700) resonance as a narrow peak 
around 50 MeV in the photon spectrum from the radia- 
tive decay of i/'(3770) into V(3700)7. 

The predicted radiative decay, obtained following an 
approach which has been successful in predictions of the 
radiative decay of into ao(980) and /o(980), should 
be rather reliable. Based on the calculated rate, the 
planned run of 0^(3770) production at BEPC-II should 
produce enough events to gather sufficient statistics of 
decays into the V (3700)7 channel. The study done here 
should stimulate experimental efforts in this direction. 
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